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Abstract. - The scaling behavior of self-avoiding walks (SAWs) on the backbone of percolation 
clusters in two, three and four dimensions is studied by Monte Carlo simulations. We apply the 
pruned-enriched Rosenbluth chain-growth method (PERM). Our numerical results bring about the 
estimates of critical exponents, governing the scaling laws of disorder averages of the end-to-end 
distance of SAW configurations. The effects of finite-size scaling are discussed as well. 
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' Introduction. — The universal configurational prop- 
erties of long, flexible polymer chains in a good solvent 
die perfectly described by the model of self-avoiding walks 
(SAWs) on a regular lattice [1]. In particular, the aver- 
age square end-to-end distance (-R^), and the number of 
oonfigurations of SAWs with N steps obey the scaling 
Iciws: 



z 



N 



(1) 



where i^sAW and 7saw are the universal critical exponents 
tliat only depend on the space dimension d, and z is a 
non- universal connectivity constant, depending also on the 
type of the lattice. The properties of SAWs on a regular 
lattice have been studied in detail both in computer sim- 
ulations [2-7] and analytical approaches [8-11]. For ex- 
ample, in the space dimension d=3 one finds within the 
frame of the field-theoretical renormalization group ap- 
proach i'sAW=0.5882 ± 0.0011 [11] and Monte Carlo sim- 
ulations give i'sAW=0.5877± 0.0006 [5]. For space dimen- 
sions d above the upper critical dimension (iup=4, the scal- 
ing exponent becomes trivial: ^'sAw(rf > 4)=l/2. 

A question of great interest is how SAWs behave on 
randomly diluted lattices, which may serve as a model of 
linear polymers in a porous medium. Numerous numer- 
ical [12-22] and analytical [17, 23-29] studies lead to the 
conclusion that weak quenched disorder, when the con- 
centration p of lattice sites allowed for the SAWs is higher 
than the percolation concentration pc, does not influence 
the scaling of SAWs. The scaling laws (1) are valid in this 



case with the same exponents, independently of p. More 
interesting is the case, when p equals the critical concen- 
tration Pc (see Table 1) and the lattice becomes percola- 
tive. Studying properties of percolative lattices, one en- 
counters two possible statistical averages. In the flrst, one 
considers only incipient percolation clusters whereas the 
other statistical ensemble includes all the clusters, which 
can be found in a percolative lattice. For the latter en- 
semble of all clusters, the SAW can start on any of the 
clusters, and for an A^-step SAW, performed on the ith 
cluster, we have (i?^) If, where li is the averaged size of 
the ith cluster. In what follows, we will be interested in 
the former case, when SAWs reside only on the percolation 
cluster. In this regime, the scaling laws (1) hold with new 
exponents ^ t^sAW,7pc 7SAW [14-16,20,21,30,31]. 
A hint to the physical understanding of this phenomenon 
is given by the fact that weak disorder does not change 
the dimension of a lattice, whereas the percolation cluster 
itself is a fractal with fractal dimension d^^ dependent on d 
(see Tabic 1). In this way, scaling law exponents of SAWs 
change with the dimension d^^ of the (fractal) lattice on 
which the walk resides. Since dup=G for percolation [39], 
the exponent Vp^{d > 6)=l/2. For the connectivity con- 
stant Zp^ of SAWs on a percolative lattice the estimate 
Zp^=PcZ is suggested, where z is the value on the corre- 
sponding pure lattice [32] . 

Until recently there did not exist any satisfactory the- 
oretical estimates for scaling law exponents of SAWs on 
percolation clusters, based on refined field-theoretical ap- 
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Table 1: Critical concentration pc of site-diluted lattices and 
fractal dimensions of percolation cluster dp^ and backbone of 
the percolation cluster d^^ for different space dimensions d. 



d 


Pc 


dp. 


< 


2 


0.592746 [33] 


91/49 [36] 


1.650 ± 0.005 [38] 


3 


0.31160 [34] 


2.51 ± 0.02 [37] 


1.86 ±0.01 [38] 


4 


0.19688 [35] 


3.05 ± 0.05 [37] 


1.95 ±0.05 [38] 



proachcs. In particular this was caused by the rather com- 
plicated diagrammatic technique of the perturbation the- 
ory calculations. Recently the field theory developed by 
Meir and Harris [17] was reconsidered in Refs. [30,31], 
where the field theory with complex interacting fields has 
been constructed and a special diagrammatic technique 
developed. The scaling properties of a SAW on a percola- 
tion cluster were found to be described by a whole family 
of correlation exponents z/^*-', with i^'-^'^z^sAW- 

Note that up to date there do also not exist many studies 
dedicated to Monte Carlo (MC) simulations of our prob- 
lem and they do still exhibit some controversies. The first 
MC study of a SAW statistics on a disordered (diluted) 
lattice in three dimensions was performed in the work of 
Kremer [12]. It indicates no change in the exponent v 
for weak dilution, but for concentrations of dilution near 
the percolation threshold a higher value i/p^ w 2/3 was 
observed. 

This result was the only numerical estimate of Vp^ for 
a number of years, until Lee et al. [13, 14] performed MC 
simulations for a SAW on the percolation cluster for square 
and cubic lattices at dilutions very close to the percola- 
tion threshold. Two earlier of these works indicate the 
rather surprising result that in two dimensions the critical 
exponent is not different compared to the pure lattice 
value. Later, some numerical uncertainties were corrected 
and the value for Vp^ found in two dimensions is in a new 
universality class. This result has been confirmed in a 
more accurate study of Grassberger [15]. In the case of 
three and four dimensions, there also exist estimates in- 
dicating a new universality class [14], but no satisfactory 
numerical values have been obtained so far. It was argued 
in Ref. [20] , that series enumerations of all possible SAW 
configurations on a percolation cluster give a greater value 
for Vp^ (and therefore in better agreement with theoretical 
prediction) than that obtained from MC simulations due 
to some specific peculiarities of the latter method. 

In the present paper, the so-called chain-growth algo- 
rithm is applied. Conventional MC methods such as mul- 
ticanonical sampling [40] or the Wang-Landau method [41] 
expose problems in tackling "hidden" conformational bar- 
riers in combination with chain update moves which usu- 
ally become inefficient at low temperatures, where many 
attempted moves are rejected due to the self-avoidance 
constraint. Rosenbluth chain growth avoids occupied 
neighbors at the expense of a bias. Chain-growth methods 



Fig. 1: Percolation cluster on a d=2-dimensional regular lattice 
of edge length L=50. 



with population control such as PERM (pruned-cnrichcd 
Rosenbluth method) [42, 43] improve the procedure con- 
siderably by utilizing the counterbalance between Rosen- 
bluth weight and Boltzmann probability. PERM has been 
applied successfully to a wide class of problems, in partic- 
ular to the 0-transition of homopolymcrs [42] , trapping of 
random walkers on absorbing lattices [44] , study of protein 
folding [45], etc. 

Construction of percolation clusters. We con- 
sider site percolation on regular lattices of edge lengths up 
to Lmax=300, 200, 50 in dimensions d=2, 3, 4, respectively. 
Each site of the lattice was assigned to be occupied with 
probability Pc (values of critical concentration in different 
dimensions are given in Table 1), and empty otherwise. 
To extract the percolation cluster, we apply the algorithm 
of site labehng, based on the one proposed by Hoshen and 
Kopelman [46]. If for a given lattice it is not possible to 
find a cluster that wraps around in all coordinate direc- 
tions, this disordered lattice is rejected and a new one is 
constructed. The typical structure of percolation clusters 
is presented in Fig. 1. On finite lattices the definition 
of spanning clusters is not unique (e.g., one could con- 
sider clusters connecting only two opposite borders), but 
all definitions are characterized by the same fractal dimen- 
sion and are thus equally legitimate. The here employed 
definition has the advantage of yielding the most isotropic 
clusters. Note also that directly at p = Pc more than one 
spanning cluster can be found in the system, and the prob- 
ability P{k) for at least k separated clusters grows with 
the space dimension as P{k) exp(— afc''/^'^"^-') [47,48]. 
In our study, we take into account only one cluster per 
each disordered lattice constructed, in order to avoid pre- 
sumable correlations of the data. 

Since we aim on investigating the scaling of SAWs on 
a percolative lattice, we are interested rather in the back- 
bone of the percolation cluster, which is defined as follows. 
Assume that each bond (or site) of the cluster is a resistor 
and that an external potential drop is applied at two ends 
of the cluster. The backbone is the subset of the cluster 
consisting of all bonds (or sites) through which the cur- 
rent flows; i.e., it is the structure left when all "dangling 
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Fig. 2: Elastic (left) and geometrical (right) backbones of the 
percolation cluster depicted in Fig. 1. 



ends" are eliminated from the cluster. The SAWs can be 
trapped in "dangling ends" , therefore infinitely long chains 
can only exist on the backbone of the cluster. 

The algorithm for extracting the backbone of a given 
percolation cluster was first introduced in [49] and im- 
proved in [50] . This so-called burning algorithm is divided 
into two parts. First, we choose the starting point -"seed" 
- at the center of the cluster. For all the sites on the edge 
of the lattice, belonging to the percolation cluster, we find 
the shortest paths between the "seed" and the given site. 
As a result, we obtain a so-called skeleton or elastic back- 
bone [51] of the percolation cluster, shown in Fig. 2, left. 

In the second part of the algorithm, wc consider succes- 
sively each site of the elastic backbone and check whether 
a "loop" starts from this site. A "loop" is a path of sites, 
belonging to the percolation cluster, which is connected 
with the elastic backbone by no less than two sites. Sites 
of the elastic backbone together with sites of "loops" form 
finally the geometric backbone of the cluster (see Fig. 2, 
right). 

Once a cluster is generated, its fractal dimension in 
topological (or chemical) space I can be determined ac- 
cording to [50]: 

{Mb{1)) - l'- , (2) 

where Mb{1) is its "mass" (number of cluster sites) and 
is the fractal dimension of the backbone in chemical 
space. It is related to the dimension d^^ in coordinate 
—dgdntin, where dmin is the fractal dimension 



space by dp^ = 

of the shortest path on the backbone and describes the 
scaling behavior between r and I, i.e. (l) ^ j-^mia^ with 
d„,i„=1.130 ± 0.004 in d=2 [52], (imin=-1.374 ± 0.003 in 
d=3 [34], (i,nin=1.567 in rf=4 [53]. The results for fractal 
dimensions of the percolation cluster and its geometrical 
backbone in d=2,3,4 are compiled in Table 1. 

The method. — We use the pruned-enriched Rosen- 
bluth method (PERM), proposed in the work of Grass- 
berger [42] . The algorithm is based on ideas from the very 
first days of Monte Carlo simulations, the Rosenbluth- 
Rosenbluth (RR) method [2] and enrichment strategies 
[54]. Let us consider the growing polymer chain, i.e., the 



nth monomer is placed at a randomly chosen neighbor site 
of the last placed {n — l)th monomer {n < N, where N 
is the total length of the chain). In order to obtain cor- 
rect statistics, if this new site is occupied, any attempt 
to place a monomer at it results in discarding the entire 
chain. This leads to an exponential "attrition" , the num- 
ber of discarded chains grows exponentially with the chain 
length, which makes the method useless for long chains. In 
the RR method, occupied neighbors are avoided without 
discarding the chain, but the bias is corrected by means 
of giving a weight Wn ^ iY\"=2 '^0^^ to each sample con- 
figuration at the Tith step, where m; is the number of free 
lattice sites to place the Ith monomer. When the chain 
of total length N is constructed, the new one starts from 
the same starting point, until the desired number of chain 
configurations are obtained. The configurational averag- 
ing for the end-to-end distance r = y/B?{N) is then given 
by: 



(r) = 



TT/conf ^conf 
conf ^^N ^ 



wconf 

Qf ^'^ N 



:^rP(r, N), 



(3) 



where 1^^°"^ is the weight of an A^-monomer chain in a 
given configuration and P{r, N) is the distribution func- 
tion for the end-to-end distance. 

While the chain grows by adding monomers, its weight 
will fiuctuate. PERM suppresses these fiuctuations by 
pruning configurations with too small weights, and by 
enriching the sample with copies of high-weight config- 
urations [42]. These copies arc made while the chain 
is growing, and continue to grow independently of each 
other. Pruning and enrichment are performed by choos- 
ing thresholds W^f and depending on the estimate 
of the partition sums of the n-monomer chain. These 
thresholds are continuously updated as the simulation pro- 
gresses. The zeroes iteration is a pure chain-growth algo- 
rithm without reweighting. After the first chain of full 
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Fig. 3: Disorder averaged distribution function rP(r, N) vs the 
scaling variable r/A^"'''^ in d=3 dimensions. Lattice size L=200, 
number of SAW steps A=40 (squares), A=50 (pluses), iV=60 
(diamonds), N=70 (crosses), A''=80 (stars). 
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Fig. 4: Disorder averaged distribution function rP{r, N) vs 
the scaling variable r/N'^^" in d—A dimensions. Lattice size 
L—50, number of SAW steps A'^=15 (squares), A'^=18 (trian- 
gles), iV=20 (pluses), N=25 (crosses), iV = 30 (stars). 



length has been obtained, we switch to W^, . If the 
current weight Wn of an n-monomer chain is less than 
, a random number r=0, 1 is chosen; if r=0, the chain 
is discarded, otherwise it is kept and its weight is dou- 
bled. Thus, low- weight chains are pruned with probabil- 
ity 1/2. If Wn exceeds , the configuration is doubled 
and the weight of each copy is taken as half the origi- 
nal weight. For updating the threshold values we ap- 
ply similar rules as in [43,45]: W^=C(Z„/Zi)(c„/ci)2 
and W^=0.2Wn , where c„ denotes the number of cre- 
ated chains having length n, and the parameter C con- 
trols the pruning-enrichment statistics. After a certain 
number of chains of total length is produced, the it- 
eration is finished and a new tour starts. We adjust the 
pruning-enrichment control parameter such that on aver- 
age 10 chains of total length A^ are generated per each 
iteration [45], and perform 10^ iterations. Also, what is 
even more important for efficiency, in almost all iterations 
at least one such a chain was created. 

Results. — To study scaling properties of SAWs on 
the backbone of percolation clusters, we choose as the 
starting point the "seed" of the cluster, and apply the 
PERM algorithm, taking into account, that a SAW can 
have its steps only on the sites belonging to the backbone 
of the percolation cluster. In the given problem, we have 
to perform two types of averaging: the first average is per- 
formed over all SAW configurations on a single backbone 
according to (3); the second average is carried out over 
different realizations of disorder, i.e. over many backbone 
configurations: 

_ 1 c 



Table 2: Results of linear fitting of obtained results for (r) 
for SAWs in d—3 dimensions on the backbone of percolation 
clusters, L=200. denotes the sum of squares of normalized 
deviation from the regression line, DF is the number of degrees 



of freedom. 


mill 




a 


X'/DF 


6 


0.665 ± 0.003 


0.946 ± 0.003 


2.783 


11 


0.668 ± 0.003 


0.935 ± 0.004 


2.269 


16 


0.669 ± 0.003 


0.930 ± 0.004 


2.054 


21 


0.669 ± 0.003 


0.924 ± 0.004 


1.345 


26 


0.667 ± 0.002 


0.930 ± 0.006 


0.743 


31 


0.668 ± 0.002 


0.934 ± 0.008 


0.844 



Here, C is the number of different clusters, the index c 
means that a given quantity is calculated on the cluster c, 
P(r, N) is the distribution function, averaged over cluster 
configurations. 

The case of so-called "quenched disorder" is considered, 
where the average over different realizations of disorder 
is taken after the configurational average has been per- 
formed. As it was pointed out in [15], the correctness of 
results, obtained in the picture of "quenched" disorder, 
depends on whether the location of the starting point of 
a SAW is fixed while the configurational averaging is per- 
formed, or not. In the latter case, one has to average over 
all locations and effectively this corresponds to the case of 
annealed disorder. Thus, as we have already stated above, 
we start each configuration of a SAW on the same site - 
the "seed" of the backbone of a given percolation clus- 
ter. We use lattices of the size up to Lniax=300, 200, 50 in 
(i=2,3,4, respectively, and performed averages over 1000 
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Fig. 5: Disorder averaged end-to-end distance vs number 
of steps in double logarithmic scale for SAWs on the back- 
bone of percolation clusters in d—2 (pluses), d=3 (stars), d=4 
(squares). Lines represent linear fitting, statistical error bars 
are of the size of symbols. 



Scaling behavior of self-avoiding walks on percolation clusters 



4 


y' — 


4.5 
4 






Table 3: The exponent v^^ for a SAW on a percolation clus- 
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3: Averaged end-to-end distance vs number of steps 
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our results 


0.782 ± 0.003 


0.667 ± 0.003 
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The disorder averaged distribution function (5) can be 
written in terms of the scaled variables r/ (r) as: 



rP(r, N) - /(r/(r)) ^ /(r/A^^^- )■ (6) 

The distribution function is normalized according to 
P{r, N)^l. The numerical results for the distribu- 
tion function in d=3 and d=4 are shown in Figs. 3 and 
4 for different N. When plotted against the scaling vari- 
able r/N^P" , the data are indeed found to nicely collapse 
onto a single curve, using our values for the exponent 
reported in Table 3 below. 

To estimate the critical exponents i^p^, linear least- 
square fits with lower cutoff for the number of steps A^min 
are used. The value (sum of squares of normalized de- 
viation from the regression line) serves as a test of the 
goodness of fit (see Fig. 5 and Table 2). 

Since we can construct lattices only of a finite size L, 
it is not possible to perform very long SAWs on it. For 
each L, the scaling (1) holds only up to some "marginal" 
number of SAWs steps A^marg, as it is shown in Fig. 6. We 
take this into account when analyzing the data obtained; 
for each lattice size we are interested only in values of 
A^ < A'niarg, which results in effects of finite-size scaling 
for critical exponents. 
Let us assume that N, 



marg 



L'^ , and for a SAW confined 



inside a lattice with size L finite-size scaling holds: 



N 



where 



g{x)— const when N 



(7) 



so that Eq. (1) is recovered. The crossover occurs at 
(r) ~ L, N=N^s_rg, which leads to uj~1/i^. Here, v stands 
for vsAW or i/p^ for the cases of the pure lattice and back- 
bone of percolation cluster, respectively. Similar scaling 
properties have already been observed in problems of ran- 
dom walks in confined environment in Ref. [55]. 



Having estimated values for the critical exponent I'p^, 
presented in Table 3, we can proceed with testing the 
finite-size scaling assumption (7). When plotted against 
the scaling variable N/L'^, the data for different lattice 
sizes L should collapse onto a single curve if we have found 
the correct values for the critical exponents. The numer- 
ical results for the scaling function g{N/L'^) both for the 
pure lattice (for comparison) and the backbone of perco- 
lation clusters are presented in Fig. 7. Note, that our 
estimation of the exponent i^sAW in two dimensions gives 
0.745 ± 0.002. 

Conclusions. — The present paper concerns the uni- 
versal configurational properties of SAWs on pcrcolative 
lattices. The statistical averaging was performed on the 
backbone of the incipient percolation cluster, which has 
a fractal structure and is characterised by fractal dimen- 
sion dp^. Note, that up to date there do not exist many 
works dedicated to Monte Carlo (MC) simulations of our 
problem and they do still exhibit some controversies. In 
particular, in the case of four dimensions, there exist only 
estimates, indicating a new universality class [14], but no 
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Fig. 7: The scaling function g[N / L^) as a function of its 
argument at data collapse for three different lattice sizes 
L=100, 150, 200 in d=2. Left: pure lattice, right: backbone 
of the percolation cluster. 
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satisfactory numerical values for critical exponents have 
been obtained so far. 

Applying the pruned-enriched Rosenbluth method 
(PERM), we studied SAWs on the backbone of percolation 
clusters, using lattices of size up to Lniax=300, 200, 50 in 
d— 2,3,4, respectively, and performing averages over 1000 
clusters in each case. Our results bring about numerical 
values of critical exponents, governing the end-to-end dis- 
tance of SAWs in a new universality class. The effects of 
finite lattice size are discussed as well. 



V.B. is grateful for support through the "Marie Curie 
International Incoming Fellowship" EU Programme and to 
the Institut fiir Theoretische Physik, Universitat Leipzig, 
for hospitality. 

REFERENCES 



DES Cloizeaux J. and Jannink G., Polymers m Solu- 
tion (Clarendon Press, Oxford) 1990; de Gennes P.-G., 
Scaling Concepts in Polymer Physics (Cornell University 
Press, Ithaca and London) 1979. 

Rosenbluth M.N. and Rosenbluth A.W., J. Chem. 
Phys., 23 (1955) 356. 

Madras N. and Sokal A.D., J. Stat. Phys., 50 (1988) 
109. 

Mac Donald D., Hunter D.L., Kelly K. and Naeem 
J., J. Phys. A, 25 (1992) 1429. 

Li B., Madras N. and Sokal A.D., J. Stat. Phys., 80 
(1995) 661. 

Caracciolo S., Causo M.S. and Pelissetto A., Phys. 
Rev. E, 57 (1998) R1215. 

MacDonald D., Joseph S., Hunter D.L., Moseley 
L.L., Jan N. and Guttmann A.J., J. Phys. A, 33 (2000) 
5973. 

J.C. and ZiNN-JuSTiN J., Phys. Rev. B, 21 



[2] 

[3] 
[4] 
[5] 
[6] 
[7] 

[s; 

[9 

[lo; 
[11 

[12 
[13 



[14 

[is: 

[16; 

[17 

[is; 

[19 

[2o; 

[21 
[22; 



Le Guillou 
(1980) 3976. 

NiENHUIS B. 

Le Guillou 



Phys. Rev. Lett, 49 (1982) 1062. 
J.C. and ZiNN- Justin J., J. Physique Lett., 
46 (1985) L127; J. Physique, 50 (1988) 1365. 
GuiDA R. and Zinn Justin J., J. Phys. A, 31 (1998) 
8104. 

Kremer K., Z. Phys. B, 45 (1981) 149. 
Lee S.B. and Nakanishi H., Phys. Rev. Lett, 61 (1988) 
2022; Lee S.B., Nakanishi H. and Kim Y., Phys. Rev. 
B, 39 (1989) 9561. 

Woo K.Y. and Lee S.B., Phys. Rev. A, 44 (1991) 999. 

Grassberger p., J. Phys. A, 26 (1993) 1023. 

Lee S.B., J. Korean Phys. Soc., 29 (1996) 1. 

Meir Y. and Harris A.B., Phys. Rev. Lett, 63 (1989) 

2819. 

Lam P.m., J. Phys. A, 23 (1990) L831. 
Nakanishi H. and Moon J., Physica A, 191 (1992) 309. 
RiNTOUL M.D., Moon J. and Nakanishi H., Phys. Rev. 
E, 49 (1994) 2790. 

Ordemann A., Porto M., Roman H.E., Havlin S. and 

BuNDE A., Phys. Rev. E, 61 (2000) 6858. 

Nakanishi H. and Lee S.B., J. Phys. A, 24 (1991) 1355. 



[23 
[24; 

[25 
[26 

[27 
[28 

[29 

[3o; 

[31 

[32 

[33 
[34' 
[35 

[36 

[37 

[38 

[39 

[40 

[41 

[42 
[43' 

[44; 

[45; 

[46 

[47 
[48; 

[49 

[50 

[51 

[52 

[53; 

[54; 

[55 



Kim Y., J. Phys. C, 16 (1983) 1345. 

Barat K., Karmakar S.N. and Chakrabarti B.K., J. 

Phys. A, 24 (1991) 851. 

Sahimi M., J. Phys. A, 17 (1984) L379. 

Rammal R., Toulouse G. and Vannimenus J., J. Phys. 

(Pans), 45 (1984) 389. 

Kim Y., J. Phys. A, 20 (1987) 1293. 

Roy A.K. and Blumen A., J. Stat Phys., 59 (1990) 

1581. 

Lam P.m. and Zhang Z.Q., Z. Phys. B, 56 (1984) 155. 
von Ferber C, Blavatska V., Folk R. and Holo- 
VATCH Yu., Phys. Rev. E, 70 (2004) 035104(R). 
Janssen H.-K. and Stenull O., Phys. Rev. E, 75 (2007) 
020801(R). 

Chakrabarti B.K., Bhadra K., Roy A.K. and Kar- 
makar S. N., Phys. Lett A, 93 (1983) 434. 
ZiFF R.M., Phys. Rev. Lett, 72 (1994) 1942. 
Grassberger P., J. Phys. A, 25 (1992) 5867. 
Paul G., Ziff R.M. and Stanley H.E., Phys. Rev. E, 
64 (2001) 026115. 

Havlin S. and Ben Abraham D., Adv. Phys., 36 (1987) 
155. 

Grassberger P., J. Phys. A, 19 (1986) 1681. 
Moukarzel C, Int. Journal of Mod. Phys. C, 8 (1998) 
887. 

Stauffer D. and Aharony A., Introduction to Percola- 
tion Theory (Taylor and Francis, London) 1992. 
Berg B.A. and Neuhaus T., Phys. Lett B, 267 (1991) 
249. 

Wang F. and Landau D.P., Phys. Rev. Lett, 86 (2001) 
2050. 

Grassberger P., Phys. Rev. E, 56 (1997) 3682. 

Hsu H.P., Mehra v., Nadler W. and Grassberger 

P., J. Chem. Phys., 118 (2007) 444. 

Mehra V. and Grassberger P., Physica D, 168 (2002) 

244. 

Bachmann M. and Janke W., Phys. Rev. Lett., 91 
(2003) 208105; J. Chem. Phys., 120 (2004) 6779. 
Hoshen J. and Kopelman R., Phys. Rev. E, 14 (1976) 
3438. 

Aizenman M., Nucl. Phys. B [FS], 485 ( 1997) 551. 
Shchur L.N. and Rostunov T., JETP Lett, 76 (2002) 
475. 

Herrmann H.J., Hong D.C. and Stanley H.E., J. 
Phys. A, 17 (1984) L261. 

Porto M., Bunde A., Havlin S. and Roman H.E., 
Phys. Rev. E, 56 (1997) 1667. 

Havlin S., Nossal R., Trus B. and Weiss G.H., J. 
Phys. A, 17 (1984) L957. 

Herrmann H.J. and Stanley H.E., J. Phys. A, 21 
(1988) L829. 

Janssen H.-K. and Stenull O., Phys. Rev. E, 61 (2000) 
4821. 

Wall F.T. and Erpenbeck J. J., J. Chem. Phys., 30 
(1959) 634. 

Reis F., J. Phys. A, 28 (1995) 6277. 



